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Abstract. We re-analyse earlier measurements of resistance R and piezoresistance 
K in RuCvbased thick-film resistors. The percolating nature of transport in these 
systems is well accounted by values of the transport exponent t larger than its universal 
value t ~ 2.0. Furthermore, we show that the RuC>2 volume fraction dependence of 
the piezoresistance data fit well with a logarithmically divergence at the percolation 
thresold. We argue that the universality breakdown and divergent piezoresistive 
response could be understood in the framework of a tunneling-percolating model 
proposed a few years ago to apply in carbon-black-polymer composites. We propose 
a new tunneling-percolating theory based on the segregated microstructure common 
to many thick-film resistors, and show that this model can in principle describe the 
observed universality breakdown and the divergent piezoresistance. 



PACS numbers: 72.20.Fr, 72.60.+g, 72.80.Ng 
1. Introduction 

The term thick-film resistors (TFRs) commonly indicates glass-conductor composites 
based on Ru0 2 (or also Bi 2 Ru 2 7 , Pb 2 Ru 2 6 , and Ir0 2 ) grains mixed and fired 
with glass powders of usually lead borosilicate type. The most successful commercial 
applications of TFRs exploit their quite large change in bulk resistivity when subjected 
to applied deformation. Given their low temperature sensitivity and high stability, TFRs 
are therefore often used in sensor applications like pressure or force sensors p], 0. 

Besides their applications, TFRs are also a quite unique type of metal-insulator 
composite due to their complex microstructure |3], f|, nonuniversal behavior of 
transport || |J, coexistence of quantum tunneling with percolating behavior |7[], marked 
microscopic elastic heterogeneities . These peculiarities, which to our best knowledge 
are only partially shared by other compounds, contribute to make the TFRs an 
interesting playground for basic research on metal-insulator composites. 

Concerning their microstructure, TFRs are often in a segregated structure regime in 
which large regions of glass of about 1/zm in size constraint the much smaller conducting 
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grains (from 10 /iin to about 300 //m) to be segregated in between the interstices of 
neighbouring glass grains. In addition, as a result of the firing process, there is strong 
evidence of diffussion of nanosized particles of metal insider the glass, leading therefore 
to a modified glassy region around the original metallic grains 0. 

The role played by these ultrasmall conducting particles in transport has not yet 
been elucidated. The observation that the original larger metallic grains are separated 
by thin glass films of only 1-2 nanometer thickness has suggested direct tunneling 
as the main mechanism of transport [flO . However, recent electric force microscopy 



measurements have revealed that the diffused metallic ultrasmall particles can contribute 
in a substantial way to the electrical connectivity of the conducting regions, suggesting 
therefore a possible active role in transport of such a modified glass [11], [12|] • 

In spite of their differences, the above scenarios both indicate or are at least 
consistent with tunneling transport (between adjacent original metallic grains or 
assisted by the small particles diffused in the glass). Furthermore, the tunneling 
hypothesis is sustained also by the high values of the piezoresistive response, i. e. the 
sensitivity of electrical resistance upon applied strain, and the variable-range-hopping 



like temperature dependence of transport at low temperatures [13 



The above picture, together with the percolating behavior observed as a function 
of the metallic volume concentration, suggestes that TFRs can be regarded as a 
sort of segregated tunneling-percolation systems: i. e. a non-homogeneous mixture 
of metal-insulating phases in which tunneling cohexists with percolation. A typical 
homogeneous counterpart is encountered in some carbon-black-polymer composites 
where the tunneling-percolating behavior takes place within a homogeneous and random 
arrangement of conducting and insulating phases [14|] . 

Among the different unusual aspects of TFRs, the often observed nonuniversality 
of transport is certainly one of the most interesting ones, both from the point of view 
of material science and basic statistical physics. According to the classical theory of 
transport in percolating systems, the resistance R of a metal-insulator composite with 



metallic volume concentration x follows a power-law behavior of the form [15, K 



R^Roix-Xc)-*, (1) 

where Ro is a prefactor, x c is the percolation critical concentration and t is the transport 
critical exponent. Random resistors network theories predict that Ro and x c depend on 
the microscopic details of the composite, such as the arrangement of the conducting 
phase within the sample and the elemental resistances connecting two neighbouring 
sites [15, 16]. Conversely, the value of the transport exponent t is quasi- universal: for 
percolation on a lattice of resistances with non-pathological distributions, t depends 
only upon lattice dimensionality D []I6" |. According to the most recent and accurate 
calculations, t ~ 2.0 for D = 3 jT7|, f[J. 

Deviations from universality have been repeatedly reported for various composites 
n], and, notably, for RuCvbased TFRs which can display values of t of about 
± 3-4 H |20| or even up to t ~ 7.0 |§. It is the subject of this paper to investigate 
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Figure 1. Sheet resistance Rq as a function of Ru02 volume fraction x for different 
grain sizes of the metallic powders. Open symbols: Ref. p0|; full symbols: Rcf.p3] 



further on the universality breakdown of TFRs and its microscopic origin by re- 
analysing earlier data on transport and piezoresistance. We show that reported metallic 
concentration dependences of the piezoresistive response for different Ru0 2 -based TFRs 
imply tunneling-distance dependence of the exponent t. We argue that a possible 
origin of such universality breakdown stems from tunneling-percolating processes with 
strong fluctuations in tunneling distances, as previously proposed by Balberg for D = 3 



homogeneous random arrangements of conducting and insulating phases pTJ| . Finally, 
we propose a model of TFRs where segregated structures and fluctuating tunneling 
distances are treated on the same level. 



2. Transport and piezoresistance of Ru02-glass composites 

Among the several published experimental results concerning the transport properties 
of TFRs, many of them report the contemporary measurement of resistance and 
piezoresistance on commercial samples for which, however, the actual compositions 
are not known with sufficient accuracy. Since we are interested on the evolution of 
transport as the conducting concentration is varies in a controlled way, the literature 
based on commercial TFRs is of little use. However, we have been able to single out two 



published sets of data of Ru0 2 -based TFRs |22| which, to the best of our knowledge, 
are the only reports providing both the sheet resistance Ra and the piezoresistance as 
a function of the Ru0 2 volume concentration x. According to Refs. f20|, p2[ , both sets 
of measurements have been done for Ru0 2 powders with different specific surface areas 



(SA) mixed with borosilicate |^tj or high lead silicate |22[] glasses, and the final resistors 
were obtained after standard firing cycles. 

In figure (0) we report the sheet resistance R\j data as a function of Ru0 2 volume 
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Figure 2. Universal log-log plot relating (R^/Rq) 1 ^ with the distance x — x c from the 
critical concentration. Open symbols: Ref . p0| ; full symbols: Ref . p2[ . The straight 
line is equation ([!]). Inset: t vs x c obtained by the fitting of equation ([j]) to the 
experimental data. The dotted lines are a guide to the eye. 



fraction for different Ru02 grain sizes $ as extracted from the SA values ||23|| . The 
open symbols refer to the data reported in Ref. |20|| , while the full symbols to Ref. ||22|| . 
All sets of data show a strong increase of Rq as x decreases, reflecting the vicinity 
to a percolation thresold where Rn diverges. The percolating nature of transport in 
these composites is demonstrated in the log-log plot of figure (H) where all data collapse 
in a single straight line in accord with the power-law behavior of equation (|l|). The 
values of Rq, the critical concentration x c and of the exponent t have been obtained by 
least-square fits to (H), and are reported in table [l|. 

There are a few points worth to be discussed. First, as it is clear from the data 
of Ref. [20] in table [I], the critical concentration x c decreases as the surface areas SA 
(or mean grain size $) of the RuC>2 powders increases (decreases). This trend is also 
observed in the data of Ref. |22 and it can be understood by realizing that finer RuC>2 
powders occupy more likely the interstitial regions between neighbouring glass grains. 
As shown by Refs.|| f|, [|, this enhanced segregation effect is reflected in a reduction 
of x c . Note that, however, the SA vs x c data of Ref. p0| do not merge with those of 
Ref. [^2] . Maybe the different nature of the glass used in these two types of TFRs is at 
the origin of this discrepancy. 

A second important point is that the values of the transport exponent t are well 
above the universal value t ~ 2.0 of three dimensional lattices. This result confirms 
other similar findings |3], || and clearly indicates breakdown of transport universality in 
RuC>2-based TFRs. It is worth to remind that transport nonuniversality implies that 
the exponent t acquires a dependence on microscopic details such as local intergrain 
resistances and microstructure. In this respect, we note that the samples of Refs. |20|, [22] 
display an increase of t with the critical volume fraction x c as shown in the inset of figure 
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Table 1. Specific surface areas (SA) and corresponding RuC>2 grain sizes $ for the 
samples used in Refs. ||(], . Ro, x c , and t are the values of the parameters of equation 
(Q) extracted by least square fits of the data reported in Refs.j2(J [2^]. Kq and B values 
are the corresponding fits of the piezoresistance data to equation (^). The bracketed 
values are only indicative since they have been obtained by fitting a set of only two 
points. 



SA [m 2 /g] 


$ [nm] 


Ro [fi] 


x c 


t 


K 


B 


67 a 


f3 


0.145 


0.0413 


2.71 ±0.07 


-13 ±4 


3.8 ±0.7 


13 a 


66 


0.06 


0.0592 


2.80 ±0.01 


-11 ±3 


3.8 ±0.6 


5 a 


172 


0.0392 


0.1035 


3.6 ±0.1 


-5 ±2 


3.4 ±0.5 


60 b 


14 


0.515 


0.0151 


2.42 ±0.03 


(-4.12) 


(1.73) 


9 b 


96 


0.597 


0.0581 


3.4246 


-12 ±1 


7.3 ±0.4 



a Ref. 
b Ref. 



m 

[22] 



We have noticed a qualitatively similar behavior also in the transport measurements 



of Ref. [24], where t increases from t = 2.49 to t = 5.38 varied from x c = 0.0237 

to x c = 0.0401. Whether the increase of t with general feature of Ru0 2 -glass 

compounds is still an open question and further studies are needed to settle this point. 
However it should be stressed that other ruthenate TFRs (Bi 2 Ru 2 7 and Pb 2 Ru 2 6 ) 
display a quite different behaviour, i.e., t remains constant and very close to its universal 
value (Bi 2 Ru 2 0y) or slightly decreases as x c increases (Pb 2 Ru 2 0e) [ffijfl . 

The authors of Refs. |2(| |22| have also measured the Ru0 2 concentration dependence 
of the longitudinal piezoresistance coefficient defined as: 

K L = d -M?>. (2) 
as 

The above quantity is also often called the (longitudinal) gauge factor and it is obtained 
by recording the change of the sheet resistance caused by a strain e applied along 
the direction of the external voltage drop. When e is orthogonal to the voltage drop, 
the resistance change gives rise to the transverse gauge factor K^. Both and 
depend on the geometry of the sample, but they can be expressed in terms of the 
intrinsic longitudinal and transverse piezoresistive coefficients, Ty and T±, defined as 
the logarithmic change of the resistivity under uniaxial strains applied parallel and 
orthogonal to the voltage drop, respectively P5| . The values of K L as a function of x 
are reported in figure (|]) for the same samples of figures @) and (fj) f27j. Clearly, K L 
seems to diverge as x approaches to the same critical volume fraction values for which 
i?n, figure (0), goes to infinity. The origin of such a divergence has been previously 
discussed by the authors of Ref.[Q. Here, we disprove their arguments and propose an 
alternative interpretation based on the universality breakdown displayed by the data of 
figure (H) and table [|. 

Our arguments go as follows. First, since for all the x values considered the samples 
are well within the critical region \x — x c \ 1, then there is not marked difference 
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Figure 3. Values of the longitudinal piezoresistance -Kl as a function of RuC>2 volume 
fraction x for different specific areas of the metallic powders. Open symbols: Rcf . p(i| ; 
full symbols: Rcf.||. 



between longitudinal and transverse piezoresistive responses, and the samples can be 
considered, in a first approximation, as electrically isotropic even if an uniaxial strain 
has been imposed to evaluate K^. This is due to the tortuosity the current path has in 
flowing through the sample when x is close to x c . More specifically, the piezoresistive 
anisotropy factor defined as x = (T|| — T±)/T\\ goes to zero as (x — x c ) x with A ~ 0.5 for 
three dimensional systems | 26fl , so that for \x— x c \ <C 1 and apart from geometric factors, 
the x-dependence of uniaxial strain is basically indistinguishable from that obtained by 
applying equal strain along all directions. Hence, in the critical region, it is a good 
approximation to evaluate K L directly by differentiation of equation ([!]) with respect to 
e, just as done in Ref . [pCf| : 

Kia ^_^ Hx xc)] (3) 

The first term is the usual contribution to the piezoresistance and it generally 
measures the strain sensitivity of the mean junction resistance between two neighbouring 
conducting particles. Since, by construction, this term is independent of the Ru0 2 
concentration, all the x dependence of Ki, reported in figure (^) must come from the 
last term of equation (P) . The authors of Ref. |2D|| have assumed that t and x c are not 
strain sensitive, and argued instead that x, being a volume concentration, depends on e. 
In this way, from equation (|3]), they obtained a (x — a; c ) _1 divergence of for x —>■ x 
which fitted in a satisfactory way their data [open symbols in figure (§)]. There is 
however a misconception in this reasoning. In fact, in percolation theory, x is actually a 
measure of the concentration p of (tunneling) intergrain junctions with finite resistances 
present in the sample. Current can flow from one end to another of the composite as 
long as a macroscopic cluster of junctions spans the entire sample. Hence the application 
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Figure 4. Universal semi- log plot of (Xi — Kq)/B. The straight line is equation 



of a strain, which is typically only of order e ~ 10 -4 , can eventually change the value of 
the junction resistances (affecting therefore Ro) but cannot modify the concentration p 
of junctions |28| . 



The dependence of x upon strain claimed in Ref . [g0| is therefore erroneous and the 
x dependence of Kj, must have a different origin. Since also x c is e- independent, the only 
quantity left which could have a strain dependence is the exponent t. This hypothesis is 
perfectly in agreement with the universality breakdown demonstrated by the t data of 
table 0. Hence, by allowing a nonzero derivative of t with respect to e we obtain from 
equation (|3|): 

1 



K L ~ K + B In 

where we have defined 

d\n(R ) 



x 



K 



de 



B 



dt 
de' 



(4) 



(5) 



In figure we report the results of fitting the data of figure to equation (f|), 
represented by the straight line. The agreement is fairly good and indicates that the 
data are indeed consistent with a logarithmic divergence at x — * x c . The fitted values 
of the parameters Kq and B are reported in the last two columns of table |I] where 
the main feature is that B is always positive while Kq is always negative. This latter 
result is quite interesting. In fact, naively, one expects that if Rq is a measure of the 
mean interparticle tunnel junction resistance, then it is a good approximation to set 
Ro qc exp(2a/£), where a is the mean tunneling distance and £ oc 1/DV is the tunneling 
factor (or localization length) and V is the tunneling barrier potential. Under a tensile 
strain e, the mean tunneling distance is enhanced by ae, while £ remains constant. 
In this situation therefore K ~ 2a/£ is always positive, contrary to the K values of 
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table |I[ We shall see in the next section that this apparent contradiction is solved by 
allowing for sufficiently strong fluctuations in the tunneling distances. In this case in 
fact the notion of mean tunneling distance a can be meaningless and Kq can assume 
also negative values. 



3. Balberg's model 

In the previous section, the logarithmic divergence of has been obtained by allowing 
a strain-dependence of the resistance exponent t. From a microscopic point of view, 
this property is equivalent to require that t is a functional of the intergrain tunneling 
distances r tun between two neighbouring Ru0 2 particles. In fact, an applied strain e 
would modify the intergrain distance r tun to r tU n(l +e) leading therefore to a variation of 
t. Actually, a microscopic theory of this kind exists and has been proposed by Balberg 
a few years ago as a possible mechanism of nonuniversality in carbon-black-polymer 



composites ||21|| . Briefly, the theory of Balberg goes as follows. Consider a random 
resistor network where the interparticle (bond) conductance distribution is given by: 

p(g)=ph(g) + (l-p)5(g), (6) 

where p is the fraction of bonds with finite conductance g with distribution h(g). The 
resistor network is constructed in order to mimic a system constituted by a random 
arrangement of spherical particles of given diameter $ and it is assumed that electrons 
can tunnel from one particle to its nearest neighbouring particle with a tunneling 
junction conductance of the form: 

g = g e~ 2r ^, (7) 

where go is a constant, £ is the localization length assumed to be smaller than the average 
interparticle distance a, and r tun = r — $ is the distance between the surfaces of two 
neighbouring spheres whose centers are separated by r. Since the particle centers are 
distributed more or less randomly, the distance r depends on the particular couple of 
particles considered. Hence to obtain the distribution function h(g) of the tunneling 
conductances one needs the distribution function P{r) of inter-particles distances. 
Balberg proposed that the relevant distribution function is of the form: 

p -r/a 

P{r) = — , (8) 

Lb 

where a is the mean inter-sphere distance and a>$. By using 

POO 

%)= / drP(r)5(g-g e- 2r ^), (9) 
then the distribution function of tunneling conductances reduces to 

m = l=2.(S-Y, (10) 



'1 1 ^ 



go \go 

where 

1 £ 
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is supposed to vary between a = and a = 1. As it is well known, distributions of the 
form of equation fllPf ) lead to nonuniversal behavior of the transport exponent t ||29| . In 
fact, according to renormalization group results and highly disordered resistor networks 
analysis, the overall resistance close to the percolation thresold is given by equation ([I]) 
with IBtt HI: 



t 



(12) 



t Q if {d-2)v +!/{!- a) < t 

(d - 2)v + 1/(1 - a) if (d - 2)v + 1/(1 - a) > t 

where to is the universal value of the exponent (t — 2.0 in three dimensions), d is the 
dimension of the system, and v is the correlation-length exponent [y ~ 0.88 in three 
dimensions). We have therefore arrived at the result that, as soon as a is sufficiently 
large, the transport exponent depends upon the tunneling distance a and so can be 
affected by an applied external strain e. Hence, from equations ([]]), (|[), and (|i~2l) it is 
easily found that the piezoresistance K L follows equation (Q) with 

B = £~- (13) 
de 1 — a 

Since a < 1, the parameter B is always positive, in agreement therefore with the sign 
of the B values reported in table |l[ 

The distribution function of equation fllO|) can also explain the negative values 
of Kq shown in the last column of table [l| and that, as we have pointed out in the 
last section, are not expected when ordinary interparticle conductance distributions are 
considered. In fact, let us consider an effective medium approximation (EM A) to the 
bond conductance distribution of equations @ and ([T0|): 

G-g 



dgp(g)- 



o. 



g + 2G 

where G is the effective macroscopic conductance [TE\. For \p — p c \ <C 1 (where p. 

in EMA) the integral in equation ( p!4[ ) can be evaluated by setting x 

The resulting resistance R = 1/G is then: 

o ■ i \ 1 -l/(l-a) 
6 sm(7ra ' 



(14) 
= 1/3 

g/G with G < 1 



R 



(P-Pc 



_2 1 - a vr(l - a) 

from which the EMA piezoresistance K reduces to K 
where: 



(15) 



K + 1/(1 - a) ln(p - Pc 



K 



1 — a 



In 



3 sin(7ra) 
7r(l — a) 



+ 



7T 



1 — a tan(7ra;) 



(16) 



As shown in figure (|5|), where the above expression is plotted as a function of a, K Q is 
always negative in agreement therefore with the sign of the K values extracted by the 
fits to the experimental data of Refs. p0| , [22|j (table [1]). The possibility of having K < 
stems from the fact that the distribution function (^) has variance equal to the mean 
intergrain distance a. Hence the resistance prefactor Rq is no longer a measure of the 
mean intergrain tunnel junction resistance so that the arguments of the last section do 
not apply. 
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Figure 5. Kq in the effective medium approximation [equation (|T6|)] as a function of 
the tunneling parameter a = 1 — £/2a. 



The above results for B and Kq have been derived by a straightforward derivative 
of i? with respect to the applied strain e. However, elastic heterogeneities can have 
important effects. This should be particularly true for Ru0 2 -glass systems where quite 
stiff metallic particles are embedded in a much softer glass. The main effect of such 
elastic heterogeneity is that the external applied strain induces highly varying local 
strains which depend on the relative bulk moduli of the metallic and insulating phases 
and on the microstructure ||. We can approximately describe the main effect of elastic 
heterogeneity by arguing that the tunneling distances are affected by a local strain 
E\ oc ~ Ae, where A is a function of the relative bulk moduli of the conducting and 
insulating phases. In this way, B and K of equation @ reduce to: 

K = A^^, B = A-p—. (17) 

Since the local strains are concentrated mainly within the softer (insulating) phase 
through which the electrons tunnel, then A is expected to be larger than the unity ||. 
The elastic heterogeneity has therefore the main effect of amplifying the piezoresistive 
response by enhancing the absolute values of B and Kq, however leaving their signs 
unchanged. 



4. Tunneling and percolation in segregated compounds 

In the previous sections we have provided evidences that RuCVbased thick-film resistors 
have nonuniversal behavior of transport and that their piezoresistance response fits well 
with a logarithmic divergence about the critical concentration x c . Furthermore, we 
have shown that such a divergent piezoresistance arises naturally from a specific model 
of tunneling transport in percolative composites proposed some years ago by Balberg. 
Although these results are quite encouraging, we are going to see that there are several 
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Figure 6. tunneling conductance distribution function of equation (19) resulting from 
a random arrangement of conducting spheres in three dimensional space. 



open questions regarding Balberg's theory in general and its applicability to TFRs in 
particular. 

Let us start by considering the interparticle distribution function of equation @. 
The reasons which prompted Balberg to use equation @ are discussed in his original 
where the interested reader can directly refer. 



paper 



21 



However, if we literally 
consider a three dimensional random arrangement of spheres, then instead of equation 
(0) the inter-sphere distance distribution function for low volume densities is more 



correctly given by the Hertz distribution |32 
P(r 



3r(4/3) 3 ^e' 



-r(4/3) 3 (r/a)3 



where T(4/3) ~ 0.893. As shown in Ref.[[33|, the above expression is also the exact 
distribution function for a random set of fully impenetrable spheres at arbitrary density. 
For large r, equation ([18]) falls to zero more rapidly than equation (§), leading to a quite 
drastic effect on h(g). In fact, from equations (0), (H) and and again assuming for 
simplicity $ -C a, the resulting tunneling conductance distribution function reduces to: 



Kg) 



3r(4/3) 3 (l 



a 



In 2 



l-r(4/3) 3 (l-a) 3 In 2 (g /g) 



(19) 



9o\ 9_ 
9o \9 ) V£o, 

The above expression is plotted in figure |6] for different values of a. Clearly, in 
contrast to equation (|T0"D, h(g) does not diverges for g — > but instead it goes to 
zero at g — 0. In Ref.f33H very accurate analytical expressions for interparticle distance 
distribution functions of impenetrable (hard-core) spheres are reported. We checked 
that h(g = 0) = also for this case so that such a vanishing limit is a general 
feature of randomly distributed spheres in three dimensional space. Since transport 
nonuniversality and the appearance of a tunneling distance dependence of the exponent 
t have as necessary condition a divergent h(g) behavior for g — > 0, then the distribution 
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function of conductances resulting from a random arrangement of (impenetrable or 
not) spheres is not able to account for nonuniversality and the consequent logarithmic 
divergence of the piezoresistive response. Such a conclusion is quite at odds with 
Balberg's theory and its extension for the piezoresistive problem discussed in the 
previous section. However there are a couple of points worth to be stressed. First, 
as shown in figure [], h(g) has a maximum which shifts at lower values of g and sharpens 
as a is made larger. Already for a = 0.6, the peak of h(g) occurs at g/g ~ 0.03 <C 1. In 
this situation, the critical region where the resistance follows equation ([!]) with universal 
exponent t = t — 2.0 is expected to be considerably narrowed. Hence if p is not very 
close to p c the resistance could appear to have the form of equation ([[]) with an apparent 
nonuniversal value of t. 

A second important point concerns the use of interparticle distance distribution 
functions derived from random arrangements of spheres. In real composites, interactions 
between conducting and insulating phases and microstructural inhomogeneities can 
easily make the idealized model of randomly distributed spheres of little use. A 
remarkable example is provided actually by thick-film resistors where the segregation 
effect due to the different sizes of glass and RuC>2 particles is completely missed by, for 
example, equation fll8|). However the segregation effect has interesting consequences on 
the tunneling distance distributions. Let us consider the microstructure model of TFRs 
proposed originally by Pike || (see also Refs.||, In this model, the conducting 
particles (approximated by spheres of diameter $) are arranged in such a way to occupy 
narrow channels along the edges of large cubes of side L>$ which represent the glass 
particles. These channels form a cubic lattice spanning the whole sample and each 
channels can be occupied or free of conducting spheres. We can imagine the channels 
to be quasi-one dimensional so that each occupied channel is actually given by a series 
of adjacent metallic particles separated by the sintered glass. Each channel can contain 
up to about L/<3> number of spheres and, for typical TFRs, L/$ is roughly between 
50 and 500 depending on the relative sizes of the glass and conducting particles. The 
interesting point here is that equation (§) is the exact interparticle distance distribution 
function for a one-dimensional random arrangement of penetrable spheres. Hence, h(g) 
given by equation fllOj) is the relevant conductance distribution function for the junction 
conductances inside an occupied channel. The question is now whether the total channel 
conductance G has a power law diverging distribution as well. If a channel is occupied 
by L/$ > N ^> 1 spheres, then the total channel conductance is that given by N 
junctions in series each with conductance (i — 1, . . . , N): 

N 1 1 - 
G- 1 = E - - 9 m l + ~(N - 1)0X (20) 

where the distribution function flip] ) has been used and g m i n is the minimum interparticle 
conductance for a row of N particles. In writing equation (p0|), we have set go = 1 and 
have implicitly assumed that N is large enough to replace the summation with an 
integral. The quantity g m i n (and so G) varies from channel to channel and to find its 
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Figure 7. Computer generated distribution functions F(G) of the total conductances 
G of a series of N tunnel junctions with conductances distributed according to equation 
@ with g = 1 and a = 0.5. The dotted curves are F(g) oc G~ Q b . 



distribution function f(g m i n ) consider the probability that, in a given channel, g min is 
larger than a given arbitrary value 5: 

Prob(# min > 6) = / dg min f(g min ). (21) 



The above probability must equal the probability that none of the N junctions within 
a channel has conductance smaller than 5: 

,5 nJV 



Prob(# min > 5) 



1 - / dgh(g) 



(l_ S i-<*} N . (22) 



/(<?min) is then readily found by equating the right hand sides of equations (|2lf) and (^2]) 
and by taking the derivative with respect to 5: 

/(ft-ta) = (1 - a)tf(l " ^^"te (23) 
The final distribution function of channel conductances G is then obtained by 

combining equation ( |2"0D with equation (PB|). An asymptotic estimate of F(G) for small 
G can be derived from equation (ED) by realizing that for g min (l/^) 1-0 the total 
channel conductance is dominated by the smallest junction conductance: G ~ (?mm- 
Hence in this limit and from equation (p3|): 

F(G?) oc G a , for G < (l/iV) 1 -". (24) 

The validity of the above limit is demonstrated in figure |^ where we have plotted 
computer generated values of F(G) for a = 0.5 and for different number iV of junctions. 
Clearly, in accord with equation (Ejj), the numerical distribution function diverges as 
G~ a independently of N. 

The set of channels which form the cubic lattice in our model have therefore 
conductances G with diverging distribution at small G. In this situation therefore, 
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by following the arguments of the previous section, transport universality breaks down 
and logarithmically divergent piezoresistive response sets in as long as a is sufficiently 
large. This interesting result can be made even more stringent by considering the more 
realistic situation in which the spheres inside the channels are impenetrable. Following 



Ref. [33 1, the exact interparticle distance distribution function for this case is: 



P (r) = _L^ e -(r-*)/(«-*) j (25) 

where r > <3> and, again, $ and a (a > $) are the diameter of the spheres and the mean 
interparticle distance, respectively. By using equations (|7|) and (f|) and by following the 
above analysis, it is straightforward to find that the distribution function h(g) within 
a channel and the channel conductance distribution F(G) are again given by equations 



([i"0|) and (|24"D, respectively, where however a of equation (|TT|) is now replaced by: 



& = l--&-. (26) 
a — <P 

This result has not required the hypothesis that r ^> $, as in the original formulation 
of Balberg, and has therefore a more general validity. A more complete analysis of 
our model including Monte Carlo calculations and a generalization to spheres with 
distributed diameters will be presented elsewhere. 

Before concluding this section, it is worth to mention that in modelling the 
microstructure of TFRs along the lines of Pike [H] we have implicitly considered the 
metallic particles as more or less the original grains used in fabricating the TFRs. We 
have therefore not addressed the role of the small metallic clusters of nanometer size 
which have been shown to occupy the space between two neighbouring RuC>2 grains 
as result of metal-glass interactions |9|, [U], [12]. It could be argued that these particles 
effectively reduce the tunneling potential barrier separating the metallic larger grains, 
enhancing therefore £. Further investigations and possible modifications of our model 
are therefore necessary for a more complete description of transport in TFRs. 

5. Conclusions 

In summary, we have confirmed that the sheet resistance of Ru02-based thick-film 
resistors |20], |22| follows equation (|I|) with transport exponent t larger than the universal 



value t ~ 2.0 characteristic of three dimensional systems [17], |i8| . Furthermore, we have 
shown that the piezoresistance data taken on the same samples [20, 22] are consistent 
with a logarithmic divergence close to the percolation thresold. We have argued that 
such divergence stems from a tunneling-distance dependence of the transport exponent t 



and an analysis of Balberg's theory |TJ, of tunneling-induced universality breakdown 
seems to qualitatively reproduce some of the experimental findings of Refs. [p0|, [22| . 
Furthermore, we have adapted Balberg's theory, which was originally proposed for 
homogeneously random arrangements of conducting particles in an insulating host, to 
the highly segregated structure of TFRs. We have shown that this model naturally gives 
rise to a diverging conductance distribution function leading to universality breakdown 
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of transport and logarithmically divergent piezoresistive response at the percolation 
thresold. 
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